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We study a new mechanism for the production of dark matter in the universe which does not rely 
on thermal equilibrium. Dark matter is populated from the thermal bath subsequent to inflationary 
reheating via a massive mediator whose mass is above the reheating scale, Trh- To this end, we 
consider models with an extra U(l) gauge symmetry broken at some intermediate scale (Mint — 
10 10_12 GeV). We show that not only does the model allow for gauge coupling unification (at a 
higher scale associated with grand unification) but can naturally provide a dark matter candidate 
which is a Standard Model singlet but charged under the extra U(l). The intermediate scale gauge 
boson(s) which are predicted in several E6/SO(10) constructions can be a natural mediator between 
dark matter and the thermal bath. We show that the dark matter abundance, while never having 
achieved thermal equilibrium, is fixed shortly after the reheating epoch by the relation T RH /Mf nt . 
As a consequence, we show that the unification of gauge couplings which determines Mi nt also fixes 
the reheating temperature, which can be as high as Trh — 10 11 GeV. 
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I. INTRODUCTION 

The Standard Model (SM) of particle physics is now 
more than ever motivated by the recent discovery of the 
Higgs boson at both the ATLAS Q and CMS @] detec- 
tors. The SM, however, contains many free parameters, 
and the gauge couplings do not unify. Among the most 
elegant approaches to understand some of these parame- 
ters is the idea of a grand unified theory (GUT) in which 
the three gauge couplings 0^2,3 originate from a single 
gauge coupling associated to a grand unified gauge group 
[3J. This idea is supported by the fact that quantum 
numbers of quarks and leptons in the SM nicely fill rep- 
resentations of a GUT symmetry, e.g., the 10 and 5 of 
SU(5) or 16 of 50(10). 

Another issue concerning the SM is the lack of a candi- 
date to account for Dark Matter (DM) which consists of 
22 % of the energy density of our universe. Stable Weakly 
Interacting Massive Particles (WIMPs) are among the 
most popular candidates for DM. In most models, such as 
popular supersymmetric extensions of the SM 0] , the an- 
nihilation of WIMPs in thermal equilibrium in the early 
universe determined the relic abundance of DM. 

In this letter, we will show that GUT gauge groups 
such as Eg or SO(10) which contain additional U(l) gauge 
subgroups and are broken at an intermediate scale, can 
easily lead to gauge coupling unification Q and may con- 
tain a new dark matter candidate which is charged under 
the extra U(l). However, unlike the standard equilib- 



rium annihilation process, or complimentary process of 
freeze- in Q, we propose an alternative mechanism for 
producing dark matter through interactions which are 
mediated by the heavy gauge bosons associated with the 
extra U(l). While being produced from the thermal bath, 
these dark matter particles never reach equilibrium. We 
will refer to dark matter produced with this mechanism 
as Non-Equilibrium Thermal Dark Matter or NETDM. 
The final relic abundance of NETDM is obtained shortly 
after the inflationary reheating epoch. This mechanism 
is fundamentally different from other non-thermal DM 
production mechanisms in the literature (to our knowl- 
edge). Assuming that none of the dark matter particles 
are directly produced by the decays of the inflaton during 
reheating, we compute the production of dark matter and 
relate the inflationary reheat temperature to the choice 
of the gauge group and the intermediate scale needed 
for gauge coupling unification. As an added benefit, the 
model naturally possesses the capability of producing a 
baryon asymmetry through leptogenesis, although that 
lies beyond of the scope of this work. 

The letter is organized as follows. After a summary of 
the unified models under consideration in section II, we 
show how the presence of an intermediate scale allows 
for the possibility of producing a dark matter candidate 
which respects the WMAP constraint Q and apply it 
to an explicit scenario in section III. A discussion of our 
main result is found in section IV before concluding in 
section V. 
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II. UNIFICATION IN SO(10) MODELS 

The prototype of grand unification is based on the 
SU(5) gauge group. In an extension of SU(5) one can 
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introduce SU(5) singlets as potential dark matter can- 
didates. The simplest extension in which singlets are 
automatically incorporated is that of SO(10). There 
are, however, many ways to break SO(10) down to 
SU{3) x SU(2) x U{1). This may happen in multiple 
stages, but here we are mainly concerned with the break- 
ing of an additional U(l) (or SU(2)) factor at an inter- 
mediate scale Mi n t. Here, we will not go into the details 
of the breaking, but take some specific, well-known ex- 
amples when needed. Assuming gauge coupling unifica- 
tion, the GUT mass scale, Mqut, and the intermediate 
scale Mi n t can be predicted from the low-energy coupling 
constants with the use of the renormalisation group equa- 
tion, 

H— = -fco, . (1) 



Log [at ] 




6 10 14 18 

Log[10,n] 

FIG. 1. Example of the running of the SM gauge couplings for 
SO(10) ->■ SU(4) x SU(2) L x U(X)r. 



The evolution of the three running coupling constants 
ai, ot.2 and 013 from Mz to the intermediate scale Mint is 
obtained from Eq.lJTJ) using the /3-functions of the Stan- 
dard Model: 61,2,3 = (-41/10, 19/6, 7)/2tt. We note 
that the gauge coupling, gjj, associated with U'(l) is 

related at the GUT scale to g\ of U(1)y by go = \J~\9i 

and a* = gf/Air. Between Mint and Mqut (both to 
be determined) the running coupling constants are again 
obtained from Eq.(fT]), now using /3-functions associated 
with the intermediate scale gauge group, which we will 
label bi . The matching condition between the two differ- 
ent runnings at Mint can be written: 



(oPi)- 1 +h{t mt -t Z ) 



+ bi(tint ~ t GUT ) (2) 



with t^t = In Mint, tz = InMz, t GUT = In Mqjjt, a° = 
Oi(Mz) which is measured, and a = a^Mour) is the 
unified coupling constant at the GUT scale. This gives 
us a system of 3 equations, for 3 unknown parameters: 
a, Unt, tGUT- Solving the Eq.([2]), we obtain 
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where b l3 = (bi - bj)/(bi - bj). 

To be concrete, we will consider a specific example to 
derive numerical results for the case of the breaking of 
50(10): 50(10) -> SU(4) x SU{2) L x U{l) R -+ Mint 
SU{3) C x SU(2) L x U(1) Y ->m bw SU(3) c x U(l) em . 
When the intermediate symmetry is broken by a 16 
of Higgs bosons, the bi functions are given by 61,2,3 = 
(5/2, 19/6, 63/6)/27r @, where the computation was done 
at 1-loop level. For this case, we obtain Mint = 7.8 x 10 12 
GeV and M GUT = 1.3 x 10 15 GeV using (a?^) -1 ^ 
(59.47, 29.81, 8.45). The evolution of the gauge couplings 
for this example is shown in Fig. [T] 



III. HEAVY Z' AND DARK MATTER 

It has been shown in Q and Q that a stable dark 
matter candidate may arise in SO(10) models from an un- 
broken Z 2 ~ symmetry. If the dark matter is a fermion 
(scalar) it should belong to a 3(B — L) even (odd) rep- 
resentation of SO(10). For example, the 126 or 144 
contains a stable component \ which is neutral under 
the SM, yet charged under the extra U(l). As we have 
seen, to explain the unification of the gauge couplings 
in SO(10) one needs an intermediate scale Mi nt of or- 
der 10 10 GeV. The dark matter candidate, x, can be 
produced in the early Universe through s-channel Z' ex- 
change: SM SM -> Z' -> x X- Since M z > = ^g D M mt , 
the exchanged particle is so heavy (above the reheating 
scale, as we show below) that the DM production rate is 
very slow, and we can neglect the self annihilation pro- 
cess in the Boltzmann equation. Thus while the dark 
matter is produced from the thermal bath, we have a 
non-equilibrium production mechanism for dark matter, 
hence NETDM. 



The evolution of the yield of x, Y x = n x /s follows 



dx 




{(TV) 

m x M P ^Y c 



(4) 



where n x is the number density of x and s the entropy 
of the universe, g pi g s arc the effective degrees of freedom 
for energy density and entropy, respectively; x = m x /T, 
m x being the dark matter mass, Mp the Planck mass 
and 
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X 

.(5) 

Here n eq is the equilibrium number density of the initial 
state (SM) particles; and K 1 is the first order modified 
Bessel function and k the effective degrees of freedom of 
incoming particles. 



3 



Since the production of DM occurs mainly at Trh 3> 
m x , we can neglect m x in estimating the amplitude for 
production. In this case, assuming that both x an d the 
initial state, /, are fermions, we obtain 



Mf,) 2 



(a 



s 2 (l + cos 2 



(6) 



where 9 is the angle between the two outgoing DM par- 
ticles, N[ is number of colors of the particle /, and qi is 
the charge of the particle i under U'(l) with a gauge cou- 
pling go- Here, q is an effective coupling which will ulti- 
mately depend on the specific intermediate gauge group 
chosen. With the approximations m x ,m/ -C \/s and 
Mz> 3> Tflff, and after integration over 9 and sum over 
all incoming SM fermions in the thermal bath, we obtain 
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FIG. 2. Reheating temperature as function of the SO(10) breaking 
scale for different mass of dark matter : 10, 100 and 1000 GeV 



Solving Eq. ([7J between the reheating temperature and 
a temperature T gives 

3/2 

Y X (T) 
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(8) 

the mass of the Z' by Mz' = 
-y^gDMi nt and made the approximation g p = g s . We 

note that the effect of Z' decay on the abundance of x is 
completely negligible due to its Boltzmann suppression 
in the Universe: the Z' is largely decoupled from the 
thermal bath already at the time of reheating. 

We note several interesting features from Eq. (JS)) . First 
of all, the number density of the dark matter does not de- 
pend at all on the strength of the U'(l) coupling go but 
rather on the intermediate scale (that is determined by 
requiring gauge coupling unification as we demonstrated 
in the previous section). Second, the production of dark 
matter is mainly achieved at reheating. Thirdly, once 
the relic abundance is obtained, the number density per 
comoving frame (Y) is fixed, never having reached ther- 
mal equilibrium with the bath. And finally, upon apply- 
ing the WMAP determination for the DM abundance, 
we obtain a tight constraint on Trh once the pattern of 
5O(10) breaking is known (and thus Mi nt fixed). 

Thus, given a scheme of 5O(10) breaking we can deter- 
mine the reheating temperature very precisely from the 
relic abundance constraint in the Universe. From 
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where H is the Hubble parameter and the index "0" cor- 
responds to present-day values. Combining Eq.® and 
Eq.© we find 
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TABLE I. Possible breaking schemes of SO(10). 





SO (10) -> Gx [Higgs] 


M mt (GeV) 


Trh (GeV) 


A 


4 x 2 L x 1 R [16] 


10 12.9 


3 x 10 9 


A 


4 x 2 L x 1 R [126] 


10 11.8 


1 x 10 8 


B 


4 x 2 L x 2 R [16] 


10 14.4 


3 x 10 11 


B 


4 x 2 L x 2 R [126] 


10 13.8 


5 x 10 10 


C 


3 C x 2 L x 2 R x 1 B -l [16] 


10 10.6 


3 x 10 6 


c 


3c x 2 L x 2 R x l B -L [126] 


10 8.6 


6 x 10 3 



or 



T 



RH 



2xl0 8 GeV 



Qh 2 



1/3 



/ l00GeV \ 1/3 
\ m x ) 



4/3 



10 12 GeV 7 

(11) 
= 1. We 



where we took for illustration q 2 . K 2 q 2 
show in Fig.© the evolution of Trh as function of Mi nt 
for different values of the dark matter mass m x . We can 
thus determine the reheating temperature predicted by 
different symmetry breaking pattern^]. We summarize 
them in Table HI where the values of Trh are given for 
m x = lOOGcV. 

Finally, we must specify the identity of the NETDM 
candidate in the context described above. The DM can 
be in the 126 or 144 representations of SO(10). There 
are several mechanisms to render the DM mass light , 
one of which is through a fine-tuning of the SO(10) cou- 
plings contributing with different Clebsh-Gordan coeffi- 
cients (see for example, [l(| and [ll|) to the masses of 
the various 126 components. For example, for the group 



1 We note that the value obtained for the intermediate scale in 
different 50(10) breaking schemes is not modified by the pres- 
ence of a dark matter particle which is not charged under the 
SM gauge group. 
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126(M + 2/4545^ + 2/210210^)126 (12) 

where M ~ Mgut, and a Qa singlet in a linear combina- 
tion of 21Qh and 45/r has a vev at the GUT scale. m x 
is then given by a linear combination of M and the vev 
and can be tuned to small values, while all other particles 
inside the 126 live close to Mgut- 



IV. DISCUSSION 

Unfortunately, the chance of detection (direct or indi- 
rect) of NETDM with a massive mediator Z' is nearly 
hopeless. Indeed, the diagram for the direct detection 
process, measuring the elastic scattering off a nucleus, 
proceeds through the t— channel exchange of the Z' bo- 
son, and is proportional to 1/Af|, yielding a negligible 
cross-section. In addition, due to the present low veloc- 
ity of dark matter in our galaxy (~ 200 km/s), the indi- 
rect detection prospects from s— channel Z' annihilation 
XX ~ > Z' — >• // proportional to s 2 /iVf|, is also negligible. 

As we have seen in EqjSJ the production of dark mat- 
ter occurs in the very early Universe at the epoch of re- 
heating. A similar mechanism (though fundamentally 
completely different) where a dark matter candidate is 
produced close to the reheating time is the case of the 
gravitino [U, fl2j . Indeed, in both cases equilibrium is 
never reached and the relic abundance is produced from 
the thermal background to attain the decoupling value 
T/H, with H the Hubble constant and V = (av)rif 
the production rate. However, in the case of SO(10), 
the cross section decreases with the temperature like 
(av)z' oc T 2 /M^,, whereas in the case of the gravitino 
the cross section is constant (<jv) 3 / 2 oc 1/Mp implying 
Y{T)<xT RH . 

Finally, we note that cases B and C (in Table I) predict 
reheating temperatures which are larger (B) or smaller 
(C) than the case under consideration. Case A would also 
be compatible with successful thermal leptogenesis with 



a zero initial state abundance of right-handed neutrino 
[HI, EH- However in the cases B and C, the persistence of 
the SU(2)/e symmetry would imply that the cancelation 
in Eq. [T2l would leave behind a light SU(2),r triplet (for 
DM inside a 126) or doublet (for DM inside a 144). 
These would affect somewhat the beta functions for the 
RGE's but more importantly leave behind a test of the 
model. In the triplet (doublet) case, we would expect 
three (at least two) nearly degenerate states: one with 
with charge 0, being the DM candidate, and also states 
with electric charge ±1 and ±2 (or ±1 in the doublet 
case). 

CONCLUSION 

In this work, we have shown that it is possible to pro- 
duce dark matter through non equilibrium thermal pro- 
cesses in the context of 5*0(10) models which respect the 
WMAP constraints. Insisting on gauge coupling unifica- 
tion, we have demonstrated that there exists a tight link 
between the reheating temperature and the scheme of the 
5O(10) breaking to the SM gauge group. Interestingly, 
the numerical values we obtained arc quite high and very 
compatible with inflationary and lcptogcncsis-like mod- 
els. 
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